Quasi-particle spectra of perovskites: Enhanced Coulomb correlations at surfaces 
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Photoemission spectra of the perovskites Ca^Sri-^VOa, CarcLai-^VOs, and SrRuOs indicate 
that Coulomb correlations are more pronounced at the surface than in the bulk. To investigate this 
effect we use the dynamical mean field theory combined with the Quantum Monte Carlo technique 
and evaluate the multi-orbital self-energy. These systems exhibit different degrees of band filling 
and range from metallic to insulating. The key input in the calculations is the layer dependent local 
density of states which we obtain from a tight-binding approach for semi-infinite cubic systems. 
As a result of the planar character of the perovskite t2 g bands near the Fermi level, the reduced 
coordination number of surface atoms gives rise to a significant narrowing of the surface density 
of those subbands which hybridize preferentially in planes normal to the surface. Although the 
total band width coincides with the one in the bulk, the effective band narrowing at the surface 
leads to stronger correlation features in the quasi-particle spectra. In particular, the weight of the 
quasi-particle peak near Ef is reduced and the amplitude of the lower and upper Hubbard bands 
is enhanced, in agreement with experiments. 



I. INTRODUCTION 

Strongly correlated materials are presently an active 
area of experimental and theoretical investigation 0, ■ 
Angle-resolved photoemission spectroscopy in principle 
provides a complete map of the energy and momentum 
dependent quasi-particle states and has therefore been 
used to characterize the electronic properties of many fas- 
cinating systems. As a consequence of the finite escape 
depth of the emitted electron, however, photoemission 
spectra inevitably contain information on the electronic 
structure in the bulk and near the surface. Both the 
single- and many-particle features may depend on the 
distance from the surface. For example, even in the ab- 
sence of any surface reconstruction the reduced coordina- 
tion of surface atoms causes a characteristic oscillatory 
variation of the local density of states as a function of 
the layer index, with an appreciable effective band nar- 
rowing right at the surface In strongly correlated 
materials it is precisely the width and shape of this local 
density of states which, for a given set of on-site Coulomb 
and exchange energies, determine the details of the quasi- 
particle spectrum. Accordingly, the relative weights of 
the quasi-particle peak near the Fermi level and of the 
Hubbard bands may vary with distance from the surface. 
Additional effects can arise due to reconstruction of the 
lattice at the surface, due to more pronounced rotations 
or distortions of some of the atomic groups within the 
unit cell, surface phonons, and due to enhanced electron- 
electron interaction caused by less efficient surface screen- 
ing processes. 

Surface effects in photoemission from transition metal 
oxides have been observed in several systems. Fujioka et 
al. 3 studied SrRuOs and noticed characteristic spectral 
variations caused by the frequency dependent mean free 
path of the photoelectron. A similar trend was found by 
Maiti et al. [|| for the series Ca^La^^VC^. The latter 
data suggested the highly interesting situation of a metal- 
lic bulk coexisting with an insulating surface layer. Re- 



cently, Maiti et al. and Sekiyama et al. performed 
photoemission measurements on Ca^Sri-^VOa using a 
wide range of photon energies. These data also reveal 
striking effects associated with the varying amounts of 
bulk and surface contributions to the spectra. Typically, 
the valence bands in these systems consist of a coher- 
ent peak near Ep derived from the partially filled transi- 
tion metal t2 g bands, and a satellite feature correspond- 
ing to the lower Hubbard band. Inverse photoemission 
spectra reveal an analogous Hubbard band above Ep. 
In the measurements cited above, the weight of the co- 
herent peak diminishes for shorter escape depth, while 
the satellite features (the so-called incoherent peaks) be- 
come more intense. Photoemission spectra exhibiting a 
relatively larger surface contribution therefore are more 
strongly correlated than bulk spectra. 

This explains the puzzling behavior seen in early pho- 
toemission work on Ca^Sri-^VOs H,|U E| in which Ca 
doping caused a significant suppression of intensity near 
Ep and therefore appeared to drive the system close to 
a Mott transition. These results were at odds with the 
metallic behavior found in various thermodynamic mea- 
surements independently of Ca concentration Sep- 
arating bulk and surface contributions by using differ- 
ent photon energies, Maiti et al. |fj( and Sekiyama et 
al. yj demonstrated that the bulk emission from SrVC>3 
and CaV03 is quite similar, in agreement with the low- 
frequency bulk probes and recent theoretical work |12| . 
The surface spectra of both materials, however, are con- 
siderably more correlated. 

A related example is Si'2Ru04 for which previous pho- 
toemission spectra seemed to contradict bulk de Haas- 
van Alphen measurements |13| . Recent experimental and 
theoretical work proved, however, that this discrepancy 
can be resolved by taking into account the lattice re- 
construction at the surface of Sr 2 Ru04 which leads to 
significant changes in the photoemission spectra |l4j . 

To derive reliable information on bulk properties 
of strongly correlated systems using photoemission it 
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clearly is desirable to identify single- and many-electron 
effects associated with the surface. In the present work 
we study these effects for three perovskite materials with 
widely different band fillings: the metallic compounds 
Ca x Sri_ K V03 (d 1 ) and SrRuC>3 (d 4 ), and the series 
Ca x Lai_ x V03 which is insulating for x — (d 2 ), but 
metallic for x = 0.5 (d 1 - 5 ). We evaluate the quasi-particle 
self-energy using the dynamical mean field theory based 
on the multi-orbital Monte Carlo method HE Eg. The 
important input in these many-body calculations is the 
layer dependent local density of states which we derive 
from a tight-binding scheme for semi-infinite systems Q . 
We show that the surface leads to an effective narrowing 
of the density of states of those bands hybridizing mainly 
in atomic planes normal to the surface. As a result, corre- 
lation effects at the surface are more pronounced than in 
the bulk. Such a trend had first been predicted by Pot- 
thoff and Nolting E3 who studied the metal-insulator 
phase diagram for a semi-infinite simple cubic s band at 
half filling. Here we calculate the self-energy for multi- 
band systems using realistic local densities of states for 
several cubic perovskite materials and find qualitative 
agreement with photoemission data 0, H, 0, [^Prelimi- 
nary results on SrVC-3 were published earlier |l8j . 

This paper is organized as follows. In Section II we fo- 
cus on the single-particle electronic properties of SrVC>3 
which is can be taken as representative of perovskite ma- 
terials. In particular, we discuss the evaluation of the 
layer dependent local density of states for semi-infinite 
SrV03. Section III provides the main elements of the cal- 
culation of the multi-orbital self-energy in the bulk and at 
the surface. The quasi-particle spectra of Ca2,Si'i_ x V03, 
SrRuOs, and Caa;Lai_ x V03 are presented in Section IV. 
Section V contains the summary. 



II. ELECTRONIC STRUCTURE: SRVQ 3 



In this section we discuss the bulk and surface elec- 
tronic properties of SrV03. This system can be consid- 
ered as a prototype of a cubic perovskite material since 
its one-electron structure is relatively simple, with one d 
electron per transition metal ion. The other systems can 
be qualitatively understood in terms of these properties 
by accounting for different occupations: d for SrRuOs, 
and d 2 ~ x for Ca^Lai^^VOs. 

Self-consistent electronic structure calculations for 
bulk SrVC>3 within the local density approximation 
(LDA) E3 show that the conduction bands near the 
Fermi level consist of three degenerate ti g bands derived 
from V 4+ (3d 1 ) ions. The filled O 2p bands are sepa- 
rated from the tig levels by a gap of about I eV, and the 
cubic crystal field of the V-0 octahedron shifts the V e g 
bands above the ti g bands. Because of the cubic symme- 
try, the tig bands can be represented via a tight-binding 



1.5 
1 

0.5 


-0.5 
-1 











(a) 










xz /~ 








/ yz 







M 



> 

CD 



CD 
CC 



CD 

o 




0.5 
Energy (eV) 

FIG. 1: (a) Tight-binding fit to LDA t 2g bulk bands of SrVOs 
(3d 1 ). Solid lines: dispersion within the orbital plane; dashed 
lines: dispersion perpendicular to the orbital plane (Ef = 0) . 
(b) Solid curve: density of states of SrVOs ti g bulk bands. 
Dashed curve: analogous two-dimensional density of states 
obtained by neglecting inter-planar hopping integrals. 



Hamiltonian with diagonal elements 

h xy , X y(k) = e d + t (c x + c y ) + tic x c y 

+ [*2+* 3 ( }°z , (1) 

where Cj = 2 cos(fc.;et), i = X,y, z and a is the lattice con- 
stant. Cyclic permutations yield h xz ^ xz (k) and h yz ^ yz (k). 
The ti denote effective hopping integrals representing the 
V-O-V hybridization, where to, 2, £i,3, and ti specify the 
interaction between first, second and third neighbors, re- 
spectively. For symmetry reasons off-diagonal elements 
arise only between second and third nearest neighbors 
and are of the form h xy ^ xz (k) = —t'iS y s z — t' 2 c x s y s z , 
where Sj = 2sin(fcia), i.e., they vanish at the high- 
symmetry points. Since the coefficients t[ 2 are very small 
we neglect these off-diagonal elements so that the energy 
bands are given by ej(fc) = hi t i(k), with i = xy,xz,yz. 
The tight-binding parameters and to ■ ■ ■ ti can easily 
be found by fitting the LDA energies at high-symmetry 
points of the bulk Brillouin Zone. 

Fig. I (a) shows the tig bulk bands of SrV03 along 
the main symmetry directions. The conduction bands 
in the cubic perovskite structure consist of three nearly 
non-hybridizing ti g bands. Each of these bands is ap- 
proximately two-dimensional representing weakly cou- 
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pled atomic planes. Typically, the inter-planar band 
width is about 20 times smaller than the intra-planar 
band width. To indicate the pronounced planar charac- 
ter of these bands we denote via solid lines the strong 
dispersion within the plane spanned by an orbital and 
by dashed lines the much weaker out-of-plane disper- 
sion. According to this nearly two-dimensional electronic 
structure, the bulk density of states Pb(^) exhibits the 
characteristic main peak related to the van Hove singu- 
larity at the X point of the Brillouin Zone. This is shown 
in Fig. 1 (b) where the bulk density of states is compared 
with the two-dimensional density obtained by setting the 
inter-planar hopping integrals £2,3,4 equal to zero. The 
asymmetric shape of both distributions follows from the 
second-neighbor hopping terms ~ t%. The overall shape 
of the bulk density a gree s well with the one obtained from 
LAPW calculations 

Note that as a result of the 3d 1 configuration, the 
Fermi surface of SrVC>3 consists of three nearly perfect 
intersecting cylinders containing the d xy , d xz and d yz 
states. This peculiar shape was recently observed also 
in de Haas-van Alphen measurements of CaVC>3 20] in 
spite of pronounced orthorhombic distortions. 

In the cubic environment, the three t 2g bands have 
identical density of states. At the surface, the bulk de- 
generacy is lifted since only the d xy band exhibits strong 
dispersion within the plane of the surface (the z direc- 
tion specifies the surface normal) whereas the d xz and 
dy Z bands disperse primarily within atomic planes per- 
pendicular to the surface plane. Thus, the local density 
of states of the d xy band in the first layer is similar to the 
bulk density, while that of the d xz _ yz bands is modified by 
the reduced coordination number in the z direction. To 
evaluate these surface densities we use a Green's function 
formalism Q for semi-infinite tight-binding systems. 

Neglecting again the weak hybridization between t 2g 
states, the local density of states of band i = xy,xz,yz 
is given by 

Pi, n {u) = - bnGi in (k\\,w) (2) 
fen 

where n > 1 denotes the layer index and the Green's func- 
tion can be conveniently determined from the expression 

where m = 2n - 2, \i = (4T 2 - ft) 1 / 2 and ft = w - W t . 
The parameter A denotes a surface potential. W% and Tj 
represent intra- and inter-planar contributions to the t 2g 
band energies: 

W xy = e d + t (c x + Cy) + t\c x c y 

W xz = e d + t Q c x + t 2 c y + t 3 c x c y (4) 

W yz = e d + t c y + t 2 c x + t 3 c x c y 
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FIG. 2: Layer dependent local density of states pi : „(w) of out- 
of-plane d X z : y Z bands for first three atomic planes of SrVC>3 
(Ef = 0). The surface potential A is set equal to zero. Solid 
curve: isotropic bulk density of states pb(<+>). The local density 
of the intra-planar d xy states (not shown) is similar to the bulk 
density, even in the surface layer. 

and 

T X y = t 2 + t 3 (C x + Cy) + t 4 C X Cy 

T XZ = t + tlCx + t 3 Cy + tiC X Cy (5) 

Ty Z = t + t\Cy + t 3 C x + tiC X Cy . 

How SrV03 is terminated at the surface, especially the 
charge state of the V ions in the first layer, is not yet 
known. Thus, we assume all tight-binding parameters to 
coincide with those in the bulk. The surface potential is 
chosen to ensure charge neutrality (see below) . 

Fig. 2 compares the layer dependent local density of 
states p i n {uj) of the d XZtVZ bands with the isotropic bulk 
density pb(uj). The d xz ^ yz density in the first plane is seen 
to be less sharply peaked and more narrow than Pb(uj) 
although their total widths are identical. In particular, 
Pxz,i{u) rises almost linearly below Ep in contrast to 
the plateau-like shape of pb{u). Thus, the surface spec- 
tral weight is reduced at low and high frequencies and 
enhanced at intermediate frequencies. The local density 
of d XZi y Z states in the deeper lying layers approaches the 
bulk density rather quickly, the main effect consisting in 
an oscillatory distortion of the spectral shape rather than 
any appreciable band narrowing. The rapid convergence 
towards the bulk density is to be expected because of the 
tight-binding character of the t 2g bands resulting from 
the short range of the hopping integrals. The local den- 
sity of the in-plane d xy states depends only weakly on 
the layer index and nearly coincides with the bulk den- 
sity even in the first layer (see below). 

It would of course be desirable to perform self- 
consistent electronic structure calculations for semi- 
infinite SrV03 since they should provide a more accurate 
description of the density of states in the surface region. 
Nevertheless, we believe that the key effect discussed here 
within the simplified tight-binding approach, namely, the 
preferential band narrowing of the d xz . yz states, will hold 
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quite generally. 

The single particle properties of the remaining per- 
ovskite systems discussed in this work are closely related 
to those of SrVC>3. They differ essentially by the degree 
of filling of the ti g bands, and by the deviation from cu- 
bic symmetry via orthorhombic distortions of the oxygen 
octahedra surrounding the transition metal ions. These 
distortions primarily broaden the van Hove singularity of 
the density of states but have only a minor effect on the 
overall width of the t2 g bands [l^, |21j . We discuss these 
differences in more detail in Section IV. 



III. QUASI-PARTICLE SPECTRA 

To analyze the experimental photoemission data we 
evaluate the quasi-particle spectra by taking into ac- 
count local Coulomb interactions. According to the semi- 
infinite one-electron properties discussed above we are 
dealing with a non-isotropic system where two narrow 
d X z,yz bands interact with a wider d xy band. In the 
first atomic layer this difference in effective band width is 
most pronounced and it diminishes rapidly towards the 
interior of the system. This situation is reminiscent of 
the one in the layer perovskite Sr2Ru04, which essen- 
tially consists of Ru sheets containing two nearly one- 
dimensional d XZ} y Z bands interacting with a wide intra- 
planar d xy band. The peculiar interesting feature of the 
latter system is the fact the on-site Coulomb energy lies 
between the single-particle widths of the non-degenerate 
t 2g bands: W xz , yz < U < W xy j^. In the present 
case, on the other hand, the difference between the d xz . yz 
and d xy bands at the surface is less pronounced so that 
Wi < U for all three bands. Nevertheless, since in the 
first layer the effective width of d XZjyz states is reduced, 
the effect of Coulomb correlations on the surface bands 
should be stronger than on the wider bulk bands. 

The key quantity characterizing the effect of Coulomb 
correlations on quasi-particle spectra is the self-energy 
which we evaluate here using the dynamical mean field 
theory (DMFT) 0, EE EBl A full description of cor- 
relations near the surface would be exceedingly compli- 
cated since, in principle, it would require a mixed mo- 
mentum/real space approach in order to handle the loss 
of translational symmetry normal to the surface. This 
could be accomplished using a cluster formalism in which 
the semi-infinite system is represented via a slab of finite 
thickness. Unfortunately, the planar character of the ti g 
states leads to a slowly convergent local density of d xz , yz 
states, with many spurious van Hove singularities 

stemming from the quasi-one-dimensional hopping along 
atomic rows parallel to the surface. On the other hand, a 
cluster generalization of the DMFT is feasible today only 
for very small cluster size. 

To achieve a qualitative understanding of correlations 
at SrVC>3 surfaces we ignore the momentum variation of 
the self-energy and assume that, for a given layer, it de- 
pends only on the local density of states within that layer 



[r?| . Since we neglect the weak hybridization between ti g 
orbitals, the self-energy is diagonal in orbital space. To 
evaluate the self-energy elements Ei(o>) we use the dy- 
namical mean field approach, in which is a functional 



of the bath Green's function Q i 1 
local Gi is given by 



Gi(iLUn) 



duj 



+ /i - S,(iw n ) - u)' 



where the 



(6) 



The Matsubara frequencies are denoted by ui n and fi is 
the chemical potential. On-site Coulomb correlations 
are treated using the self-consistent multiband Quan- 
tum Monte Carlo (QMC) method (for a review, see Ref. 
0)- The temperature of the simulation was 125 meV 
((3 = 8). Several runs using 64 imaginary time slices and 
I0 5 Monte Carlo sweeps were carried out. The quasi- 
particle density of states Ni(co) = — Im Gi(u)/TT was ob- 
tained via maximum entropy reconstruction [23j . 

In principle, we are now faced with a set of coupled im- 
purity problems, where the layer-dependent baths deter- 
mine the self-energies of each atomic plane. This prob- 
lem could be solved iteratively until self-consistency is 
achieved and a common chemical potential is found for 
the bulk and at the surface. Unfortunately, for multi- 
band systems this iteration procedure would be computa- 
tionally extremely demanding and will be neglected here. 
On the other hand, previous work showed that the 
self-energy is mainly governed by local Coulomb correla- 
tions and that further changes due to interlayer coupling 
are rather small. Our assumption of charge neutrality 
within the surface layer, and the fact that correlations in 
the bulk and at the surface do not differ substantially, 
most likely also reduce the importance of interlayer ef- 
fects. The key layer dependent input in our QMC-DMFT 
calculation is therefore the one-electron local density of 
states. For a given set of local Coulomb and exchange 
energies the width and shape of this local density then 
determines the details of the quasi-particle spectrum. 



IV. RESULTS AND DISCUSSION 
A. SrV0 3 

The layer dependent local density of states of SrV03 
shown in Fig. 1(b) is calculated in the absence of any sur- 
face potential. Thus, the occupation numbers also vary 
with layer index. However, since local Coulomb correla- 
tions are very sensitive to the degree of band filling, it 
seems appropriate to enforce charge neutrality by adjust- 
ing A accordingly. Fig. 3 shows the resulting local density 
of states for the surface layer. Although the main peak 
now coincides more with the bulk van Hove singularity, 
the effective narrowing is similar to the one in Fig. I (b). 
Note that the surface potential is sufficiently weak so that 
no surface states are split off below the band. Thus, the 
total band width is the same as in the bulk. No surface 
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FIG. 3: Local density of d XZtyz states of first layer of SrVC>3 
(d 1 ) in the presence of a weak surface potential A to ensure 
charge neutrality (dashed curve). The t2 g bulk density (solid 
curve) is shown for comparison. The density of the in-plane 
d xy states in the first layer (dotted curve) is nearly identical 
to the bulk density. 



potential is needed for the d xy states since their local 
density is practically identical to the bulk ti g density. 

Fig. 4 (a) shows the bulk quasi-particle density of states 
of SrVC>3 for two Coulomb energies in the region where 
the Hubbard bands which are seen as satellites in pho- 
toemission begin to emerge: U = 4.0 eV and 4.3 eV. The 
exchange energy is J — 0.7 eV |25|]. These results show 
that in the bulk U must be larger than 4 eV to obtain the 
satellite observed in photoemission spectra [|J Q ■ The 
peak near 2 eV above Ep agrees with inverse photoemis- 
sion data 01 • Although for U = 4eV there is consider- 
able correlation-induced band narrowing and an emerg- 
ing satellite shoulder, the larger U to yields an even nar- 
rower coherent feature near Ef, with the missing weight 
shifted to the lower and upper Hubbard bands. 

Note that N b (E F ) = p b (E F ) which follows (at T = 0) 
from the local approximation implicit in the DMFT for 
isotropic systems |26(. The most recent photoemission 
data [3 confirm this result. The bulk spectra shown in 
Fig. 4 (a) qualitatively agree with DMFT-LDA results for 
SrVOs by Nekrasov et al. 0] who employed somewhat 
larger values for U and J. They are also consistent with 
pre vious spectra for the 3d 1 perovskite La^Sri-^ Ti0 3 
|27| which exhibits a similar t 2g bulk density of states. 

The surface quasi-particle spectra for SrVOs are shown 
in Fig. 4(b) and (c). The lower Hubbard peak of the 
d X z,yz states in (c) is clearly visible already for U — 4eV 
because of the narrower local density of states in the 
first layer. A larger U shifts the satellite to higher bind- 
ing energies. The comparison with the spectra shown in 
Fig. 4 (a) demonstrates that correlation effects for a fixed 
value of U are stronger at the surface than in the bulk: 
The coherent peak near Ep is narrower at the surface 
and the incoherent satellite feature is more pronounced 
than in the bulk, in agreement with experiment 0, 0, • 
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FIG. 4: Quasi-particle density of states iVj(w) of SrV03 
(d ) derived from DMFT. (a) bulk ti g states; (b) surface d xy 
states; (c) surface d XZtyz states; (d) fictitious isotropic surface 
d X z,yz states (see text). Solid curves: U — 4.3 eV, dot-dashed 
curves: U = 4.0 eV (J = 0.7 eV). Dashed curve: bare densi- 
ties of states (see Fig. 3). 



d X z,yz- Although there is little single-electron hybridiza- 
tion between ti g bands, the local Coulomb interaction 
mixes them so that the d xy surface spectrum involves 
contributions arising from the more strongly correlated 



M xz,yz 



states. Another consequence of the anisotropic sur- 
face self-energy is the correlation-induced charge transfer 
between subbands [2^|. Here, we find that 0.06 electrons 
The surface quasi-particle density of d xy states in are shifted from the d xy to the d xz ^ yz bands. Also, the 
Fig. 4(b) is intermediate between Nb(oj) and N s (uj) for quasi-particle partial densities of states at Ep do not 
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need to coincide with the bare partial densities. The cou- 
pling between narrow and wide bands is a genuine multi- 
band effect and underlines the fact that single-particle 
bands in the presence of local Coulomb interactions can- 
not be considered independently. 

For illustrative purposes we show in Fig. 4 (d) the anal- 
ogous quasi-particlc spectra for a hypothetical isotropic 
system whose density of states is given by the surface 
Pxz,yz- Evidently, correlation effects are slightly more 
pronounced than at the surface of SrV03 where two such 
narrow bands interact with the wider d xy band. All spec- 
tra shown Fig. 4 are based on single-particle densities of 
identical total width but with progressively more orbital 
components exhibiting narrower spectral shape. Accord- 
ingly, the weight of the quasi-particle peak near Ef de- 
creases systematically and the upper and lower Hubbard 
bands become successively more intense. 

Note that the quasi-particle spectra in Fig. 4 are highly 
asymmetric and that the gap between the central peak 
and the lower satellite appears at smaller U than the up- 
per gap. This behavior is a consequence of the asymmet- 
ric bulk and surface density of states, and of the low band 
filling in this system (1/6 per spin band). Thus, despite 
the typical overall 3-peak structure which is characteris- 
tic of many highly correlated quasi-particle distributions, 
their is a large degree of spectral detail which is charac- 
teristic of this particular material. 

We also point out that the many-body reduction of 
the quasi-particle band width is very much larger than 
the surface-induced one-electron band narrowing. On the 
other hand, since the on-site Coulomb energy is not far 
from the critical value for a metal-insulator transition, 
the band narrowing substantially enhances the influence 
of correlations at the surface. 

It would be interesting to perform angle-resolved pho- 
toemission measurements to determine the correlation- 
induced band narrowing of the SrV03 tig bands. For 
instance the energy at T should be only a few tenths 
of an eV below Ep instead of leV as predicted by the 
LDA. Accordingly, the true ti g bands should be consid- 
erably flatter than the single-particle bands. Also, mea- 
surements using polarized light could help to separate 
correlation effects in the d xy and d XZjVZ bands. 



is quite similar to that from SrVC>3, the CaV03 spectra 
being only slightly more correlated. Thus, the new data 
are consistent with low- frequency bulk measurements. 

While SrV03 is a cubic perovskite with a V-O-V bond 
angle of 180°, the oxygen octahedra in CaV03 are dis- 
torted, so that the V-O-V bond angle is reduced to 162°, 
implying a weaker indirect hopping between d orbitals 
and a reduction of the d electron band width. As the 
most recent LDA calculations for CaV03 show this 
decrease in total band width amounts to only about 4 % 
(from 2.6 to 2.5 eV). A more striking consequence of 
the distortion is the considerable broadening of the van 
Hove singularity and the skewing of its spectral weight to 
lower frequencies. Similar results had been obtained in 
LDA calculations for SrRu03 and CaRu03 |2l|. where 
orthorhombic distortions also cause a narrowing of the 
tig bands and a broadening of the van Hove singular- 
ity. As the LDA-DMFT bulk calculations for CaV0 3 by 
Nekrasov et al. |12| show, these changes of the density of 
states lead to a small reduction of the quasi-particle peak 
near Ep and to slightly more intense Hubbard bands 
compared to those in SrV03, in agreement with the data. 

Considering the enhancement of Coulomb correlations 
at the SrV03 surface discussed in the previous section, 
and the fact that in the bulk CaV03 is already somewhat 
more correlated than SrV03, it is now plausible that cor- 
relation effects should be even stronger at the CaV03 
surface, as indeed observed. Since the local Coulomb 
energy is not far from the critical value for a metal- 
insulator transition, the tendency for enhanced surface 
correlations in CaV03 should in fact be stronger than 
for SrV03. This effect is entirely a consequence of the 
one-electron band narrowing of the d xz%yz subbands. Ad- 
ditional effects could be caused by surface reconstruction 
or by more pronounced tilting of oxygen octahedra at 
the CaV03 surface which would lead to weaker in-plane 
hopping for the d xy bands. For a more detailed analysis 
it would be necessary to carry out electronic structure 
calculations for the CaV03 surface in order to study or- 
thorhombic distortions and possible reconstructions. In 
addition, theoretical estimates of a possible increase of U 
at the surface due to reduced screening would be useful. 



B. CaVOa 

In striking contrast to thermodynamic measurements 
earlier photoemission data on Ca^Sri-^VOs 0, 0- HI] 
indicated considerably more pronounced correlations in 
CaV03 than in SrV03. These observations even led to 
speculations as to whether CaV03 might be close to at 
Mott transition. However, recent photoemission mea- 
surements taken over a wide range of photon energies 
0, demonstrated that this discrepancy can be resolved 
by carefully separating bulk and surface contributions to 
the spectra by using widely different photon energies. Ac- 
cording to these results, the bulk emission from CaV03 



C. SrRuO ;i 

The ruthenates with perovskite-based crystal struc- 
ture have attracted considerable attention during the 
past years because of a variety of fascinating properties. 
SrRuOs is the only known ferromagnetic metal among 
the Ad oxides, with a Curie temperature of 160 K. The 
one-electron density of states is similar to that of SrV03 
except for the d occupancy of the Ru derived ti g bands 
[2l|. In a perfect cubic system the Fermi level would 
nearly coincide with the van Hove singularity and the 
flat bands in this energy region. However, the oxygen 
octahedra in SrRuOs are tilted by about 8° which leads 
to a substantial broadening of the van Hove singularity. 
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FIG. 5: Solid curve: isotropic bulk density of states of cubic 
SrRuO-3; dashed curve: local density of d xz , y z states in first 
layer (E F = 0). 
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Nevertheless, the density of states near Ep remains large, 
giving rise to a large Stoner parameter. Doping with Ca 
causes further tilting of the oxygen octahedra. As a result 
the density of states at Ep is reduced even further and 
the magnetism is easily suppressed. The layered ruthcn- 
ate Sr2Ru04 is a paramagnetic metal and is currently 
intensively studied because of its unconventional p-wave 
superconductivity below 1.5 K |28j . 

The unusual magnetic properties of the ruthenates 
have stimulated wide-ranging investigations of the cor- 
relation effects in these materials. Transport measure- 
ments |s3 as we ll as theoretical studies [2^ on SrRuC>3 
seemed to suggest non-Fermi-liquid behavior. Neverthe- 
less, recent transport data on highly ordered samples [3]| 
are consistent with Fermi-liquid theory, but demonstrate 
a remarkable sensitivity to thermal and impurity scatter- 
ing. 

Photoemission spectra on SrRuC>3 Q| revealed marked 
spectral changes for different photon energies. At 
1254 eV emission from the Ru Ad bands exhibits pro- 
nounced weight at Ep and a maximum at about 1.5 eV 
binding energy. At 100 eV, on the other hand, the spec- 
tral weight at Ep is much smaller, and the intensity then 
grows down to 2 eV binding energy. (The region farther 
below is dominated by O 2p emission.) Thus, the rel- 
atively more important surface contribution caused by 
the shorter escape depth of the photoelectron makes the 
spectrum appear more correlated than at energies gov- 
erned by bulk emission. 

To investigate correlations in the bulk and at the sur- 
face of SrRuC>3 we have applied the approach discussed 
above to the case of d 4 occupancy of the ti g bands. For 
the case of cubic perovskite structure, i.e., neglecting the 
tilting of the O octahedra, the isotropic bulk density of 
states shown in Fig. 5 agrees qualitatively with results 
of LAPW calculations |2l|. The local density of d xz ^ yz 
states in the first layer shows again the characteristic shift 
of weight from low and high frequencies to the interme- 
diate range closer to Ep. The d xy surface density nearly 
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FIG. 6: Quasi-particle density of states Ni(u>) of SrRuOs 
(d 4 ) derived from DMFT. (a) bulk ti g states; (b) surface d xy 
states; (c) surface d XZtyz states; (d) fictitious isotropic surface 
d xz , yz states (see text). Solid curves: U = 3.0 eV, J = 0.2 eV; 
dashed curves: bare densities of states (see Fig. 5). 



coincides with the bulk t2 g density. 

Fig. 6 shows the quasi-particle spectra obtained from 
DMFT-QMC calculations for U = 3.0 eV and J = 
0.2 eV. As in the case of SrV03, there is a clear progres- 
sion of correlation effects as we go from the isotropic bulk 
t2 g spectra to the d xy and d XZjyz surface spectra. The hy- 
pothetical isotropic case based on the surface d XZiVZ den- 
sity (i.e., equating p xy with p xz , yz ), shown in Fig. 6 (d), is 
even more strongly correlated, with less spectral weight 
of the coherent feature near Ep and accordingly larger 
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upper and lower Hubbard bands. Again, all bare den- 
sities have the same total width. Thus, the differences 
between the bulk and surface quasi-particle spectra are 
entirely due to the different shapes of the densities. Note 
also that, in contrast to the results shown in Fig. 4 for 
SLVO3, Nb, s (EF) < pb,s{Ep)- This deviation is caused 
by the van Hove singularity at Ep and our use of a rather 
high temperature (about 1450 K) in the QMC calcula- 
tions. 

The results shown in Fig. 6 provide a qualitative expla- 
nation for the experimentally observed shift of spectral 
weight towards the lower Hubbard peak near 1.5 eV bind- 
ing energy upon lowering the photon energy from 1254 eV 
to 100 eV Q and thereby enhancing the surface contri- 
bution to the photoemission spectra. Of course, as in the 
case of SrV03, other surface effects, e.g., due to recon- 
struction, distortion of oxygen octahedra, etc., might also 
play a role and presumably would reduce not only the 
hopping between d XZtVZ orbitals of neighboring layers but 
also between d xy states within the surface plane. These 
modifications of the electronic structure would lead to a 
further enhancement of surface correlation effects. 
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The perovskite series Ca^Lai-^VOs is very interesting 
in the context of surface versus bulk Coulomb correla- 
tions since LaVC-3 [x = 0) is a Mott insulator (cP occu- 
pancy) |32| which readily becomes metallic upon increas- 
ing the Ca concentration to x = 0.5 (d 1 ' 5 occupancy). 
Since doping with Ca leads to little structural changes 
and nearly constant Hubbard U !33], the modifications 
of the electronic properties are almost entirely caused by 
the degree of band filling. 

To investigate these effects, Maiti et al. Q carried 
out photoemission measurements at Ca concentrations 
x = 0.0 . . . 0.5 and using XPS and VUV photon ener- 
gies to distinguish bulk and surface contributions. They 
observed not only the bulk Mott transition at about 
x = 0.2, but in addition a surface Mott transition at 
a slightly larger value of x. Thus, these data suggest 
the coexistence of a metallic bulk with an insulating sur- 
face layer. This finding is in conflict with predictions 01 
for the simple-cubic half-filled s band, but it is not clear 
at present to what extent the theoretical results are ap- 
plicable to multiband systems at arbitrary occupancies. 
Moreover, the calculations in Ref. were done for a 
linearized version of the DMFT in which low- and high- 
frequency excitations within the s band are essentially 
decoupled. 

We have calculated the quasi-particle spectra of LaV03 
and Cao.5Lao.5VO3 within multi-band QMC-DMFT. The 
bulk density of states is assumed to be similar to the one 
for cubic SrV03, except for the filling of the t 2g bands. 
As shown in Fig. 7, the shape of the d xz<yz density in the 
surface layer differs for these two systems since we employ 
slightly different surface potentials in the semi-infinite 



FIG. 7: Solid curves: isotropic bulk density of states of cubic 
(a) LaV0 3 (d 2 ) and (b) Cao.5Lao.5VO3 (d 1 ' 5 ). Dashed curves: 
local density of d XZiVZ states in the first layer (Ef = 0). 



tight-binding calculation to ensure charge neutrality at 
the surface. Qualitatively, however, the densities of both 
systems exhibit a similar shift of spectral weight from low 
and high frequencies to the intermediate energy range 
close to Ep. Thus, as a result of this effective narrowing 
of the local density of d xz _ yz states we can expect a similar 
enhancement of surface correlation effects as discussed 
above for SrV03 and SrRuOs. 

Fig. 8 shows quasi-particle spectra for the bulk and 
surface layer of LaVOs at Coulomb energies close to the 
metal-insulator transition. Since our QMC calculations 
are carried out at about 1450 K ({3 — 8), the resulting 
coherent peak near Ep obscures the gap between upper 
and lower Hubbard bands which would appear at low 
temperatures. Nevertheless, in view of the temperature 
dependence known from one-band systems [2], l34j the re- 
sults in Fig. 8 (a) and (b) suggest a bulk Mott transition 
in the range U — 4.5 . . . 5.0 eV. The surface spectra in 
Fig. 8 (c) and (d) demonstrate that correlation effects in 
the first layer are more pronounced than in the bulk, in- 
dicating that a gap might form at a slightly lower value 
of U. However, according to the photoemission exper- 
iments [5| LaVOs is insulating in the bulk and at the 
surface. Thus, we may conclude that within our model 
U should be about 5 eV. 

As mentioned above, doping with Ca makes LaV03 
metallic, leading to rather different excitation spec- 
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FIG. 8: Quasi-particle density of states Ni(u) of LaVC>3 
(d 2 ) derived from DMFT. (a) and (b): bulk t 2g states for 
U — 4.5 eV and U = 5.0 eV (solid curves; dashed curves: bare 
densities), (c) and (d): surface d xy states (solid curves) and 
d X z,y Z states (dashed curves) for U = 4.5 eV and U = 5.0 eV; 
J = 0.7 eV. 



FIG. 9: Quasi-particle Green's function as function of imag- 
inary time derived from DMFT (/3 — 8). (a) bulk t2 g states; 
(b) surface d X z,yz states of LaVOs (d 2 ) for U = 3.0 . . . 5.0 eV 
in steps of 0.5 eV (from below), (c) bulk t2 a states; (d) surface 
d xz , yz states of Cao.5Lao.5VO3 (d 15 ) for U = 3.0 . . . 6.0 eV 
in steps of 0.5 eV (from below). (J = 0.7 eV). Note that 
|Gi(/3/2)| is approximately proportional to Ni(Ep)- 



tra. This drastic change in electronic properties can 
be seen clearly in the different shapes of the imaginary- 
time Green's function for LaVC*3 and Cao.5Lao.5VO3, as 
shown in Fig. 9. (The value of G(r) at r « (3/2 is repre- 
sentative of the weight of the quasi-particle peak at Ep, 
while G{(3) gives the occupancy and G(0) + G(/3) = — 1 
0-) In the former case, already for U — 4.5 . . . 5.0 eV 
there is little weight near r w /3/2, indicating the reduc- 
tion of intensity near Ep. For x = 0.5, however, even 



increasing U to 6.0 eV hardly decreases G(/3/2). Thus, 
in the range of Coulomb energies where LaV03 is clearly 
insulating, i.e., U ~ 5.0 eV, Cao.5Lao.5VO3 shows robust 
metallic behavior in the bulk and at the surface, the sur- 
face spectra being slightly more correlated than in the 
bulk. Only for U = 6 eV the surface Green's function for 
Cao.5Lao.5VO3 begins to show a new feature indicating 
an instability caused by very strong local correlations. 
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The calculations discussed above suggest that LaVC>3 
is an insulator while Cao.5Lao.5VO3 is a metal, in agree- 
ment experiment. To address the interesting question 
of separate Mott transitions in the bulk and at the sur- 
face, i.e., the possible coexistence of metallic bulk and 
insulating surface properties, it is necessary to exam- 
ine the intermediate range of Ca doping concentrations, 
0.0 < x < 0.5, and to locate the critical region where the 
bulk has just become metallic while the surface is still 
insulating. This issue will be addressed in a subsequent 
publication. 

V. CONCLUSION 

In summary, we have performed multiband QMC- 
DMFT quasi-particle calculations for several perovskite 
materials using realistic local densities of states in the 
bulk and at the surface. As a result of the planar na- 
ture of the ti g states, there is an appreciable effective 
narrowing of the local density of states in the first layer, 
even though its total width coincides with the one in 
the bulk. This band narrowing is entirely a consequence 
of the reduced atomic coordination at the surface and 
does not depend on the existence of extra surface states. 
Typically such states are due to split-off states below or 
above the main conduction band and would not neces- 
sarily contribute to the band narrowing. Other effects, 
however, which we have not included here, such as tilt- 
ing or distortion of oxygen octahedra in the surface layer, 
might well lead to a further reduction of in-plane and/or 
out-of plane d — d hopping parameters and to additional 
band narrowing. We also note that the calculations dis- 
cussed in the present work assume the same value for the 
Coulomb and exchange energies at the surface as in the 
bulk. In fact, reduced screening processes close to the 
surface could cause an increase of U and thereby make 
surface correlations even more significant. Moreover, the 
termination of the perovskite surfaces, in particular, the 
charge state of the transition metal ions, should be in- 
vestigated. 

The main point of our work is that the effective nar- 



rowing of the surface local density of states leads to an 
enhancement of correlation effects compared to those in 
the bulk, in agreement with photoemission data on vari- 
ous systems. Since photoemission spectra always involve 
bulk as well as surface contributions, it is clearly im- 
portant to identify the latter in order not to associate 
them erroneously with enhanced bulk correlations. The 
pronounced two-dimensional characteristics of d states 
near Ep is a common feature of many transition metal 
oxides. Thus, the surface-induced enhancement of corre- 
lation effects discussed in the present work should be a 
phenomenon observable in many materials. 

We finally remark that all of the systems discussed in 
the present work have very similar, albeit highly asym- 
metric bulk densities of states. They differ primarily by 
the degree of band filling, i.e., by the position of the 
Fermi level with respect to the density of states features. 
Also, the surface densities of the out-of-plane subbands 
of the various systems are qualitatively similar. Never- 
theless, even though the quasi-particle spectra for suf- 
ficiently large on-site Coulomb energies all exhibit the 
familiar three-peak structure, with a quasi-particle peak 
at Ep and lower and upper Hubbard satellites, the po- 
sitions, widths, asymmetric shapes and relative weights 
of these three peaks vary considerably for the different 
materials. To describe these variations which evidently 
reflect the original densities and band fillings, one would 
need at least 11 parameters (the total area integrates 
to unity), suggesting that even in the strong-correlation 
limit the spectra reveal a high degree of system-specific 
information concerning their electronic structure. Obvi- 
ously, far greater complexity can be expected for non- 
cubic and magnetic materials. The adequate description 
of this complexity which lies at the root of many fas- 
cinating material properties underlines the tremendous 
progress that has been achieved by combining density 
functional theory and quantum impurity methods via the 
Dynamical Mean Field Theory. 
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